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ABSTRACT 

Airborne and space-based low-resolution spectroscopy in the 1980s discovered tantalizing 
quantitative relationships between the gas phase C/O abundance ratio in planetary nebulae 
(PNe) and the fractions of total far-infrared luminosity radiated by the 7.7 and 11.3 fim bands 
(the C=C stretch and C-H bend, respectively), of polycyclic aromatic hydrocarbons. Only a 
very small sample of nebulae was studied in this context, limited by airborne observations 
of the 7.7 fim band, or the existence of adequate IRAS Low Resolution Spectrometer data 
for the 11.3 fim band. To investigate these trends further, we have expanded the sample of 
planetaries available for this study using Infrared Space Observatory (ISO) low-resolution 
spectra secured with the Short Wavelength Spectrometer (SWS) and the Long Wavelength 
Spectrometer (LWS). The new sample of 43 PNe, of which 17 are detected in PAH emission, 
addresses the range from C/O = 0.2—13 with the objective of trying to delineate the pathways 
by which carbon dust grains might have formed in planetaries. For the 7.7-/im and 1 1.3-/im 
bands, we confirm that the ratio of band strength to total infrared luminosity is correlated with 
the nebular C/O ratio. Expressed in equivalent width terms, the cut-on C/O ratio for the 7.7- 
/Lim band is found to be O.60 4, in good accord with that found from sensitive ground-based 
measurements of the 3.3-/i band too. 
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1 INTRODUCTION 

Emission features occur near 3.3, 6.2, 7.7, 8.7, 11.3 fj,m in in- 
frared (IR) nebular spectra (Russell, Soifer & Merrill 1977; Russell, 
Soifer & Willner 1977; Sellgren, Werner & Dinerstein 1983; Aitken 
& Roche 1982; Roche & Aitken 1986). Airborne spectroscopy (Co- 
hen et al. 1986,1989) shows these form a generic spectrum, with 
the most intense peaks unobservable from the ground (6.2,7.7 /im); 
all band intensities are correlated; and the 7.7-/im band tightly 
correlates with gas phase C/O ratio in planetaries (i.e., for newly 
formed dust, in a circumstellar environment). The features are at- 
tributed to vibrationally excited polycyclic aromatic hydrocarbons 
(PAHs) and related materials (Duley & Williams 1981; Leger & 
Puget 1984; Allamandola, Tielens & Barker 1985; Peelers et al. 
2002). The PAHs are regarded as permeating almost every phase of 
the interstellar medium (ISM)(Allamandola, Hudgins, & Sandford 
1999). Various components are recognized (Allamandola, Tielens 
& Barker 1989; Tielens 1993): the narrow features are carried by 
molecular-sized PAHs (50 C-atoms); larger PAH clusters (500 C- 
atoms) yield the plateaus underlying the narrow features; while the 
25 and 60-/im "cirrus" emission originates from amorphous carbon 
grains, perhaps built from PAHs, of size 5000 to 50,000 C-atoms. 
This hypothesis is supported by the obvious link between carbon 
abundance and the bands (Barlow 1983; Cohen et al. 1986,1989; 



Casassus et al. 2001), and the tight correlation between the 6.2 and 
7.7-^m bands (both from C=C skeletal modes). The 11.3-/^m fea- 
ture is due to out-of-plane bending vibrations of peripheral H atoms 
attached to aromatic molecular units; its precise frequency depends 
on the number of adjacent H atoms on each edge ring (Bellamy 
1958). The observed 11.3 and 12.7-/im bands and the 10.5-14 ^m 
plateau of emission (Cohen, Tielens & Allamandola 1985) are char- 
acteristic of aromatic rings with nonadjacent, or 2 or 3 adjacent, 
peripheral H atoms (see also van Diedenhoven et al. 2004). With 
the advent of ISO (Kessler et al. 1996) spectra, additional features 
attributed to PAHs have been identified, such as the 16.4-/im peak 
and plateau, attributed by Moutou et al. (2000) to C-triple bonds. 

Using exclusively the IRAS Low Resolution Spectrometer 
(LRS) data base, Volk & Cohen (1990) sought a correlation of 
the 11.3-/im band with C/O, finding an apparently constant value 
of 11.3-/im strength for C/0>1, with an abrupt transition near 
C/O = 1, the boundary between O- and C-rich nebulae. Any rela- 
tionship between the 1.1-^m or 1 1.3-/^m bands and gas phase C/O 
may also hold clues to the mechanisms of dust formation in plan- 
etaries, and might even elucidate the role of peripheral hydrogen 
atoms. 

Understanding carbon dust formation has wide significance in 
astronomy, because AGB stars, the precursors of C-rich PNe, are 
the dominant formation sites known for refractory grains 
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Table 1. ISO SWS and LWS data of PNe gathered for this study 
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Table 1 - continued 



Object 


TDT 


AOT 


Date 


Integration(s) 


LWS off? 


Diameter(") 


Diameter reference 


NGC6445 


48700507 


SWSOl 


17-Mar-1997 


1912 




33 


Zhang &Kwok (1993) 


NGC6445 


48700508 


LWSOl 


17-Mar-1997 


640 


Y 






NGC6537 


47000722 


SWSOl 


28-Feb-1997 


1912 




75 


Schwarz et al. (1992) 


NGC 6537 


47000802 


LWSOl 


28-Feb-1997 


1318 


Y 






NGC 6537 


70300475 


SWSOl 


18-Oct-1997 


3454 








NGC 6543 


02400714 


SWSOl 


ll-Dec-1995 


6544 




20 


Zhang & Kwok (1993) 


NGC 6543 


02400807 


SWSOl 


ll-Dec-1995 


3484 








NGC 6543 


02400910 


SWSOl 


ll-Dec-1995 


1096 








NGC 6543 


02800908 


SWSOl 


15-Dec-1995 


1094 








NGC 6543 


03201202 


SWSOl 


19-Dec-1995 


1860 








NGC 6543 


25500701 


LWSOl 


28-Jul-1996 


1266 








NGC 6572 


30901603 


LWSOl 


21-Sep-1996 


1266 


Y 


14. 


Zhang & Kwok (1993) 


NGC 6572 


31901125 


SWSOl 


30-Sep-1996 


1140 








NGC 6720 


17601005 


LWSOl 


ll-May-1996 


4321 


N 


70 


Zhang & Kwok (1993) 


NGC 6720 


36600206 


LWSOl 


16-NOV-1996 


1268 


N 






NGC 6720 


36600207 


SWSOl 


16-NOV-1996 


1140 








NGC 6741 


13401806 


SWSOl 


30-Mar-1996 


1062 




7.8 


Zhang & Kwok (1993) 


NGC 6790 


13401107 


SWSOl 


30-Mar-1996 


1062 




1.8 


Zhang & Kwok (1993) 


NGC 6790 


13401608 


LWSOl 


30-Mar-1996 


1330 


N 






NGC 6826 


30201113 


LWSOl 


14-Sep-1996 


1266 


N 


25 


Zhang & Kwok (1993) 


NGC 6826 


30201114 


SWSOl 


14-Sep-1996 


1140 








NGC 6884 


13901709 


SWSOl 


04-Apr-1996 


1834 




6.0 


Zhang & Kwok (1993) 


NGC 6886 


13400810 


SWSOl 


30-Mar-1996 


1062 




7.4 


Zhang & Kwok (1993) 


NGC 6891 


37600943 


SWSOl 


27-NOV-1996 


1912 




12.6 


Zhang & Kwok (1993) 


NGC 6891 


37600944 


LWSOl 


27-NOV-1996 


1166 


Y 






NGC 7009 


34400517 


LWSOl 


25-Oct-1996 


1268 


N 


27 


Zhang & Kwok (1993) 


NGC 7009 


34400518 


SWSOl 


25-Oct-1996 


1140 








NGC 7009 


73801242 


SWSOl 


23-NOV-1997 


1912 








NGC 7027 


02401183 


SWSOl 


ll-Dec-1995 


6542 




15 


Acker et al. (1992) 


NGC 7027 


23001256 


SWSOl 


04-Jul-1996 


1140 








NGC 7027 


23001357 


SWSOl 


04-Jul-1996 


1912 








NGC 7027 


23001458 


SWSOl 


04-Jul-1996 


3453 








NGC 7027 


55800537 


SWSOl 


27-May-1997 


6537 








NGC 7027 


71300611 


LWSOl 


28-Oct-1997 


2156 


N 






NGC 7662 


43700427 


SWSOl 


26-Jan-1997 


1912 




26 


Zhang & Kwok (1993) 


NGC 7662 


56801325 


LWSOl 


06-Jun-1997 


1318 


N 






SwStl 


47101511 


SWSOl 


Ol-Mar-1997 


1140 




1.5 


Schwarz etal. (1992) 


Vy2-2 


32002528 


SWSOl 


02-Oct-1996 


1140 




0.6 


Zhang & Kwok (1993) 


Vy2-2 


54700310 


LWSOl 


16-May-1997 


1318 


N 







4 M. Cohen & M. J. Barlow 



subsequently injected into tlie ISM. Models for soot formation in 
flames have been proposed, based on neutral radicals, ions, PAHs, 
polyacetylenic chains, or fullerenes as intermediaries (Hucknall 
1985; Barnard & Bradley 1985; Curl & Smalley 1988; Ugarte 
1992,1995; Hecht 1986; Iglesias-Groth 2004). Each route could 
probably lead to C-soot, depending on specific physical conditions. 
Our goal is to explore these routes, using C-rich planetaries as labo- 
ratories to seek the precursors to C-dust grains. Strong nebular UV 
fields excite the PAHs and the resulting emission serves as a direct 
probe of the molecular gas content. 



2 THE SAMPLE OF PLANETARY NEBULAE 

There are two ways to investigate the dependence of PAHs on 
gas-phase C/0 in PNe: using a single well-studied band from the 
ground; or measuring the strongest bands from airborne and space- 
borne observatories. An example of the former approach is the sur- 
vey of the 3.3-yum band by Roche et al. (1996), who used common 
instrumentation for their entire sample of PNe. Previous space- 
based and airborne work on the relationships between the PAH 
bands and nebular C/0 has depended upon the use of measure- 
ments made with a variety of instruments, with different apertures, 
and at different spectral resolutions. This approach encounters the 
problem of poor sensitivity to either the 11.3-/xm band from the 
Kuiper Airborne Observatory (KAO) or to the 3.3-/im band from 
the ISO SWS. It is clearly important to cover a wide range of IR 
wavelengths with common instrumentation for this work, so that 
any trends found are not compromised. ISO afforded a unique op- 
portunity to secure the requisite measurements. 

However, such observations were clearly predicated primar- 
ily on the IR brightness of the nebulae, without specific regard to 
their C/0 ratios. Therefore, the PNe we have selected have well- 
determined C/0 values, with broad agreement among several au- 
thors and different analytical methodologies (e.g. Kingsburgh & 
Barlow 1994; Rola & Stasinska 1994; Zuckerman & Aller 1986). 
Further, we have isolated precisely those nebulae whose C/O val- 
ues bear most directly on the investigation of the trends with PAH 
band strength found in the old airborne data. 

Table 1 summarizes the set of PNe assembled for this study; 
objects are listed alphabetically. For each PN spectrum we tab- 
ulate the name, TDT (ISO's observation identification number), 
the astronomical observing template (AOT), date, integration time, 
whether an LWSOl spectrum was obtained in an "ofF' position, 
the angular diameter in arcsec, and a reference for the adopted di- 
ameter. We initially selected all PNe for which AOTs SWSOl and 
LWSOl (low-resolution spectra) were obtained by ISO. The pre- 
requisite for PN selection was the existence of a published value of 
C/O (eliminating such objects as M 2-43 with its bright PAH emis- 
sion but no C/O abundance ratio). The SWS aperture varies with 
wavelength. From 2-12 /im, the size is between 14" and 20" across 
the several bands; from 12-29 /jm, between 14" and 27"; and from 
29-45 jxm, between 20" and 33". The sample includes three large 
PNe that extend substantially outside the SWS apertures. One of 
these is even larger than the LWS apertures (NGC5189). For sev- 
eral nebulae of interest (i.e. with sizeable C/O) no LWSOl spectra 
were taken by ISO because of far-infrared (FIR) faintness (e.g. M 4- 
1 8) yet the PAHs were well-detected by the SWS . Therefore, we ex- 
tended the sample to PNe for which SWSOl and IRAS 60/100-/im 
photometry were available, but which lacked any usable LWSOl 
spectra (i.e. absent, or too noisy). 

The restriction to PNe for which SWSOl spectra are avail- 



able through the ISO Data Archive leads to a bimodal sample of 
nebulae. Both PNe with bright emission lines and very weak IR 
continua, and PNe of types known to exhibit PAH emission bands 
with strong IR continua (e.g. the [WCL] sequence whose central 
stars show Wolf-Rayct-likc emission spectra) were almost equally 
frequently targeted by the observing cormnunity. The integration 
times used for the sample varied by a factor of six, depending on 
the goals of the original observers. Consequently, we can offer no 
statement as to the completeness of our sample based upon any 
physical characteristic. 

For two PNe multiple data sets exist because those objects had 
been designated as LWS wavelength calibrators (NGC7027 (26 
spectra) and NGC6543 (94)). For these objects, a single represen- 
tative LWSOl spectrum was chosen. NGC6543 was also observed 
on five occasions in the SWSOl AOT but we coadded all five spec- 
tra using inverse- variance weighting to enhance the signal-to-noise 
ratio. 



3 THE AVAILABLE SPECTRA 

The two relevant ISO spectroscopic instruments are the SWS (de 
Graauw et al. 1996) and LWS (Clegg et al. 1996). In each case we 
sought low-resolution data from the AOTs "SWSOl" (covering the 
range from 2.38 to 45.2 /xm), and "LWSOl" (43-197 /xm). 

The SWSOl data include all the known mid-infrared PAH 
emission bands and emission plateaus (3.2-3.6, 6-9, 10.5-15, and 
16-21 /im) between 3 and 21 /im. The SWS apertures were gen- 
erally well-matched to the optical and radio diameters of the se- 
lected PNe. Typical achieved resolving powers were several hun- 
dred for AOTs taken at the fastest speed (i.e. shortest observing 
time ~ 1100s), although the data archive includes higher resolu- 
tion spectra for some of our chosen PNe that were observed at the 
slowest speeds (e.g. resolving power ~1500 for ~6500s observing 
time). 

LWSOl AOTs serve to assess the luminosity of the commonly 
found thermal emission from cool dust grains in PNe. The LWSOl 
spectra have a resolution of 0.3 /jm from 43-93 lira., and 0.6 jxm 
from 84-197 /im. Observing times with the LWS ranged from 640 
to 3400s, through an aperture with an effective diameter of 66-86", 
depending on wavelength. 

The Spitzer Space Telescope (SST) and its infrared spectrom- 
eter (IRS) are able to observe PNe. However, the great sensitivity of 
this observatory and its instruments implies that many of the known 
PNe with PAH emission are too bright for spectroscopy with the 
SST. Further, the small slits of the IRS cannot accommodate most 
of our target PNe. Consequently, the ISO Data Archive offers an op- 
portunity to revisit many of the PNe known to show PAH emission, 
with apertures that generally are sufficiently large to capture essen- 
tially all the PAH emission across these objects (Smith, Aitken & 
Roche 1989) or from their surrounding photodissociation regions 
(PDRs) (Aitken & Roche 1983). 

We visually inspected all SWSOl spectra for PNe that survived 
our selection criteria, seeking detections of the PAH features. Nebu- 
lae clearly showing the bands were investigated first. Subsequently 
we examined the remaining PNe with the goal of setting quantita- 
tive upper limits on their PAH emission. 
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Figure 1. Representative plots of the SWS region containing the chief PAH 
features for PNe in our sample. Dots show the individual observed points. 
The solid curves represent the smoothed spectra. 



Figure 2. Additional representative PNe SWS spectra. Details as in Fig- 
ure 



4 BAND INTEGRALS AND THE JOIN OF SWS AND LWS 
SPECTRA 

The PAH bands are low-resolution features. Therefore, to enhance 
the detectability of the bands with SWS spectra, we interactively 
used boxcar-smoothing (with widths of 50 or 75 points) better to 
define the overall continuum and PAH features in noisy data. Fig- 
uresQand|2|illustrate the results of this smoothing applied to eight 
PNe. The IJ-fim and 11.3-^m PAH bands were integrated after 
interactively defining a single underlying continuum in each PN 
across the 5 — 15 ^im range. The wavelengths selected to define 
these continua (by cubic splines) and the regions chosen for the 
band integrals were those described by Cohen et al. (1986,1989) 
and Volk & Cohen (1990), for which choices the original relation- 
ships between integrated 7.7 or 11.3-/im intensity and C/0 were 
found. Care was taken to remove any influence of smoothed emis- 
sion lines on these band integrals. For example, the 7.46-/im Pfa 
line was cut out of those nebulae in which it was detected, by re- 
placing it by a local smooth continuum at the base of the line in the 
original, unsmoothed, SWS spectra. Likewise, in high-excitation 
PNe, we expunged the [NeVIl] line at 7.65 fim prior to smoothing 
and integrating the emission bands. In the absence of any visual 
recognition of the PAH features, formal upper limits were set by 
one of two methods. We measured the nearest positive emission 
hump lying within the wavelength interval for definition of a PAH 
integration. In addition we computed the formal integral over the 
same wavelength range of the mean spectrum plus 3 x the standard 
deviation calculated over the interval. Differences arose between 
these two methods primarily when a spectrum had negative val- 
ues or when the splined continuum locally exceeded the measured 
spectrum. We adopted whichever approach yielded the more con- 
servative upper limit. 

To evaluate the total, integrated, observed IR energy, I(IR), 
we assembled complete SWS-l-LWS energy distributions by splic- 
ing together the two data sets for each PN in the overlap region 
from 40-45 ^m. From past experience of this procedure (e.g. Co- 
hen et al. 2002) we have determined that LWSOl spectra are better 
absolutely calibrated in this region than SWS. Moreover, if a PN 
extends slightly outside the SWS aperture, the larger LWS aperture 



will more reliably assess the total nebular continuum. Therefore, 
we visually examined the overlap and rescaled each SWS spec- 
trum upwards to achieve the best match to the LWS when neces- 
sary. This approach assumes that the PN emission lying outside the 
SWS Band-4 aperture is spectrally identical to that sampled within 
the aperture. 

Given the complete spliced energy distribution, each spectrum 
was examined overall to eliminate any remaining specious features 
at the long end of the SWS band-4 range and/or the beginning of 
the LWS detector SWl range. The resulting cleaned and spliced 
spectra were then merged. When a PN had multiple spectra these 
were combined using inverse-variance weighting, before splicing 
the resulting SWS and LWS data. 

For the ten PNe lacking adequate LWSOl spectroscopy, but 
having IRAS FIR photometry, we used 60 and 100-/im flux densi- 
ties as a substitute for absent, or an alternative to noisy, LWSOl data 
to estimate the FIR contribution to I(IR). We summed the products 
of flux density and bandwidth at 60/100 /im, using the synthetic 
contiguous bandwidths of Emerson (1988) to avoid the uncertain- 
ties involved in color correction for PNe due to emission lines. Up- 
per limits at either 60 or 100 /im were treated as actual flux den- 
sities to assess the calculated FIR component of the IR luminosity. 
(These ten PNe are identified in Figures|3|and|4|by having a large 
plus sign through their plot symbols, and by asterisks in Table 2.) 
Tests were made on several PNe for which we had LWSOl data, 
comparing I(IR) derived by the two methods. These indicated that 
using IRAS data in this fashion produced values of I(IR) within 
25 percent of the actual integrated SWS-l-LWS spectra. 



5 RESULTS 

Table 2 presents the results, listing PN name, gas-phase C/O abun- 
dance ratio, I(IR) (in W cm'^), I(7.7)/I(IR), 1(1 1.3)/I(IR), and a 
reference for the C/O value adopted. Whenever possible, we have 
selected abundance ratios derived from coUisionally excited for- 
bidden lines (CELs) rather than from optical recombination lines 
(ORLs). PNe often yield substantially different abundance ratios 
when determined by the two methods (e.g. Liu et al. 2000,2001a; 
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Figure 3. Logarithmic plot for 43 PNe of tlie ratio I(7.7)/I(IR) against neb- 
ular gas-phase C/O abundance ratio. Solid line is the regression for the 15 
detected PNe after excluding CPD-56°8032 and He 2-1 13. 
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Figure 4. Logarithmic plot for 43 PNe of the ratio 1(1 1.3)/I(IR) against 
nebular gas-phase C/O abundance ratio. Solid line is the regression for the 
12 detected PNe after excluding CPD-56°8032 and He 2-113. 

Tsamis et al. 2004). Liu et al. (2000) have argued for the existence 
of H-deficient ionized clumps in PNe, cool enough to suppress for- 
bidden CELs but produce strong ORL emission. Most of the PNe in 
Table 1 have CEL C/O ratios. When these are unavailable, we have 
used the ORL line ratio (e.g., for Vy2-2). We have never mixed 
CEL and ORL abundance determinations to create a C/O value. 

In the first studies of the multi-band PAH spectra of PNe from 
the KAO it was recognised that the ratio of I(7.7)/I(IR) ranged over 
two orders of magnitude. To avoid confining the distribution of 
the bulk of the PNe detected to the lower left comer of any fig- 
ure against C/O by use of a linear plot, the logarithmic ordinate 
was selected. We have followed this precedent. CPD-56°8032 and 
He 2-1 13 are well-measured objects with extreme values of C/O, 
yet they both emit comparable fractions of the IR flux in the PAH 
bands to those emitted by other PNe with much more modest C/O 
ratios. Consequently, we exclude them from our sample in terms of 
any regression analysis. 

Fig. in illustrates the 7.7-/im results, distinguishing between 



PNe detected in the 7.7-/im band, the subset of nebulae with [WCL] 
central stars, and the 26 objects with only upper limits. The solid 
line plotted in this figure indicates the formal regression line for all 
detected PNe with the exceptions of CPD-56°8032 and He 2-113. 
The intercept is — 2.16±0.06, the slope is 1.47±0.33, and the Pear- 
son correlation coefficient for these 15 PNe is r=Q.16. We identify 
this dashed regression line as the ISO version of the relationship 
reported by Cohen et al. (1989: their Figure 20) based on airborne 
spectroscopy. (Note that, in their figure, the point for CPD-56°8032 
was misplotted too high by a factor of 10 in ordinate. It should have 
appeared at -1.1.) 

The regression line for the six [WCL] objects (after exclud- 
ing CPD-56°8032 and He 2-113) with the 7.7-^m band is insignif- 
icantly different (at the Ict level) from the line for the 15 PNe. 
M4-18 has the largest value of I(7.7)/I(IR). Ground-based spectra 
(Aitken & Roche 1982; Rinehart et al. 2002) clearly show its 11.3- 
^m band although its 8.7-^m feature is convincingly shown only 
by the Aitken & Roche (1982) spectrum from Mauna Kea. M4- 
18's SWS spectrum is very noisy, although the PAH features are 
unquestionably detected. However, the spectral levels of the PAH 
bands in this SWS spectrum match those in the spectrum by Rine- 
hart et al. (2002) to within 10 percent, confirming our estimate of 
I(7.7)/I(IR) for this object. In the small sample of PNe studied by 
Cohen et al. (1989) NGC6302 substantively helped to define the 
trend at 7.7 ^m because of its low C/O ratio. With our enlarged 
set of PNe NGC 6302 no longer controls the existence of a trend 
in this diagram. The correlation coefficient remains r=0J6 even if 
one were to exclude NGC 6302 along with CPD-56°8032 and He 2- 
113. There is no cause to reject this PN; we simply make the point 
that a single PN does not strongly influence the regression line in 
this sample at 7.7 /im. 

One immediately sees that the [WCL] PNe dominate the plot 
by contributing almost all the large ratios of I(7.7)/I(IR). Only 
NGC 6369, among the [WCL] PNe, was not detected in PAH emis- 
sion. Roche, Aitken & Whitmore (1983) observed this PN from the 
ground at very low resolving power (~40), in a 20" beam which ac- 
commodated about half the ionized gas distribution. At best there 
are possible suggestions of emission humps near 8.7 and I1.3-/im 
emission in their spectrum. The object has only a single SWSOl 
spectrum, as opposed to several of the [WCL] nebulae which have 
higher IR surface brightness. 

All the 7J-fim non-detections have relatively short exposure 
(fast speed) SWSOl observations (1062-1912 s), suggesting that 
they were spectra with poor signal-to-noise. Almost all PNe with 
C/O ratios >0.85 (logio C/O >-0.07) and I(7.7)/I(IR) >0.4 per- 
cent (logio C/O >-2.45) were detected in this PAH feature. Seven 
PNe, lying in the lower-left quadrant of Fig.|3|below these limiting 
values, were not detected in the 7.7-/im band. Only NGC 6302 is 
detected outside the above limits of C/O and I(7.7)/I(IR). This de- 
tection was certainly aided by the use of the slowest speed SWSOl 
AOT (6528 s) because the PAHs are weak in this PN and I(IR) is 
large. 

Fig. 13] similarly presents the ll3-fim results. A single re- 
gression line is again plotted for the 12 detected PNe, excluding 
CPD-56°8032 and He 2-113. The intercept is -2.59±0.08 and slope 
2. 12±0.49. The influence of NGC 6302 on this regression line is 
again minimal: r drops from 0.71 to 0.67 if one were to exclude 
this object from the regression analysis. Clearly there is no justifi- 
cation for doing this. The PN furthest from the line, with the lowest 
I(11.3)/I)IR) of the sample, is IC418 although there is likewise no 
reason to exclude it. Several nebulae show a broad hump between 
10.6 and 12.4 /im attributed to emission by SiC grains. In these ob- 
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Figure 5. Linear plot of the EW(7.7) against nebular gas-phase C/O abun- 
dance ratio. Solid hne is the regression for the 10 detected PNe for which 
meaningful EWs can be measured after excluding CPD-56° 8032 and He 2- 
113. 



jects the 11.3-^m PAH band is a very small emission feature near 
the peak of the SiC structure, requiring careful definition of a lo- 
cal continuum to extract the PAH feature. Perhaps this accounts for 
its unusual weakness in IC418, whose SiC emission is very bright. 
Overall, I(11.3)/I(IR) rises abruptly at low C/O and maintains this 
rise at least as far as C/0 ~ 2, exactly as found by Volk & Co- 
hen (1990: their Figure 3). Only the two extreme PNe are observed 
beyond this value of C/O. 

From their analysis of the 3.3-/^m band, Roche et al. (1996) 
suggested that the equivalent width (EW) of this band argued for 
a cut-off below C/O = 0.6. For high quality near-infrared ground- 
based spectra, EWs are appropriate. For space-based spectra the 
sample is still limited by signal-to-noise. This can lead to poor es- 
timates of the local continuum required to define an EW. There- 
fore, I(7.7)/I(IR) is a far more robust measure of PAH strength 
than EW for our sample of PN spectra. However, in the interest 
of trying to determine whether our sample supports a similar con- 
clusion, we have reformulated the 7.7-/im data as EWs. But this 
lack of robustness means that we must exclude several PNe from 
the sample. These are precisely those objects in which portions of 
the overall splined continua are so noisy that negative continua re- 
sult at some wavelengths within the broad 7.7-/im feature, or else 
very large EWs (>2 iim) are formally derived. We again exclude 
CPD-56°8032 and He 2- 11 3 but must also reject Sw St 1, NGC40, 
NGC6302, Cn 1-5, and M 4-18. The final sample includes 10 PNe. 
We have carried out a careful estimation of the errors in the derived 
EWs by adding in quadrature the uncertainties in each integrated 
7.7-/im intensity with those associated with the definition of the 
underlying continuum. Fractional uncertainties of 30% have been 
assigned to the observed C/0 values. Fig. |5| presents the EWs for 
these PNe linearly against C/0. The solid line is the best fitting re- 
lationship. The formal detectable onset of 7.7-/im emission is for 
C/O = 0.56o;4i, which includes the value of 0.6 found by Roche 
et al. (1996), at which C/0 value they estimate that 3.3-^m PAH 
emission first appears. 

Comparing the behaviour of CPD-56° 8032 and He 2- 1 1 3 with 
the regressions in Figures|3|and|4|suggests that perhaps these two 
objects are so extreme because some process "saturates" so that 
further increase of C/0 no longer yields more intense PAH emis- 



sion relative to the overall dust continuum emission. For both these 
PNe, essentially all of the stellar ultraviolet (UV) photons are being 
absorbed and reradiated in the IR by the PAHs and dust particles 
(Aitken et al. 1980), which is not the case for most of the other ob- 
jects in our sample. Around these two objects, there could therefore 
be a significant stratification with radius of both the spectral qual- 
ity of the ambient UV radiation field, with harder photons being 
absorbed closer in, and of the characteristic grain temperature, so 
that strong PAH emission closer in may be offset by weak PAH and 
strong far-IR grain emission further out. Therefore, the best guides 
to any underlying correlations between band strengths and C/O ra- 
tios are those provided by our PN samples that omit CPD-56°8032 
and He 2-1 13. 

The 6.2-/im PAH band is seen in 14 PNe and their relative 
strengths, I(6.2)/I(IR), are given in Table 3. In most of these neb- 
ulae this band is measured rather poorly because high noise in the 
SWS spectra below 5-6 /im adversely affects the definition of the 
underlying continua that we subtract. There is no meaningful cor- 
relation between I(6.2)/1(1R) and C/0. 



6 CONCLUSIONS 

PAHs were not detected in two of the three largest nebulae in our 
sample, either because the circumnebular PDRs that ought to con- 
tain PAHs are not sampled by the modest sizes of the SWS aper- 
tures, or else because these lower density nebulae are optically thin 
to ionizing radiation and have no PDRs. This could explain the 
absence of detectable 7.7-/im emission bands in NGC5189 and 
NGC 6720. PAHs were detected by the SWS in NGC 6537, de- 
spite its diameter of 70". However, this object is bipolar and has 
a similarly high excitation class and morphology to NGC 6302, in 
which PAHs are also detected. This morphology implies a central 
concentration of material that is lacking in large classical PNe like 
NGC 6720. Indeed, the strong 3.3 and 3.4-/im bands seen in a 3" 
slit spectrum of NGC 6537 by Roche et al. (1996) offer confirma- 
tion that the PAH emission in at least this nebula arises in the com- 
pact core. 

Using a sample of PNe more than twice as large as that pre- 
viously available from airborne spectroscopy, we have vindicated 
the existence of a relationship between I(7.7)/I(IR) and gas phase 
C/O. The fraction of total IR luminosity emitted by the 
band is observed to be roughly linearly proportional to C/0 for 
abundance ratios up to ~3. From a sample of PNe about the same 
size as that used in the LRS study by Volk & Cohen (1990) but 
based on higher signal-to-noise SWS spectra, we have confirmed 
that I(11.3)/I(IR) rises rapidly with C/0 to at least C/0 ~2, grow- 
ing roughly as the square of the gas-phase abundance ratio. For neb- 
ulae with C/O beyond a value of 2-3 this fractional growth of the 
PAH emission bands apparently ceases. The difference between the 
logarithmic intercepts of these two relationships (-2.59 and -2.16) 
indicates that, for PNe with C/0~1, the intensity of the 11.3-/im 
band is 37 percent of the 7.7-/im band's intensity. Note that the 
use of logarithmic plots greatly reduces the dependence of this ra- 
tio of band strengths on abundance ratio because the centroid of 
our sample is C/0«1. This is in excellent agreement with the av- 
erage of 36 percent measured by Cohen et al. (1986, their Table 5) 
for the generic spectrum of astrophysical PAHs. The 1986 sample 
included only 7 PNe for which band ratios relative to 1(7.7) were 
definable and included Hll regions and reflection nebulae too. We 
have, therefore, redetermined I(6.2)/I(7.7) and 1(1 1.3)/I(7.7) purely 
from our new sample of PNe, and find a mean ratio and standard 
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Table 2. Spectrally integrated observed IR fluxes, I(IR) (in W cm-^), and the ratios I(7.7)/I(IR) and 1(1 1.3)/I(IR), for 43 PNe. Names with asterisks denote 
PNe with I(IR) based on the SWS data and IRAS FIR photometry. 



PN C/O I(IR) I(7.7)/I(IR) I(11.3)/I(IR) C/O Reference 



BD+30°3639 


1.60 


4.90E-15 


0.0275 


0.0096 


Pwa, Pottasch & Mo (1986) 


CPD-56°8032 


13.1 


7.77E-15 


0.0266 


0.0056 


De Marco, Barlow & Storey (1997) 


Cnl-5 


1.7 


4.77E-15 


0.0028 


<0.0030 


Walton, Barlow & Clegg (1993) 


Hbl2* 


0.52 


1.05E-15 


<0.0027 


<0.0099 


Hyung & Aller (1990) 


He 2-1 13 


10.4 


6.00E-15 


0.0253 


0.0056 


De Marco, Barlow & Storey (1997) 


He 2-131 


0.30 


1.60E-15 


<0.0035 


<0.0015 


Liu etal. (2001a) 


Hu 1-2* 


1.23 


8.60E-17 


<0.0094 


<0.0089 


Liu et al. (2004) 


Hu2-1* 


0.48 


1.80E-16 


0.0037 


<0.0076 


Wesson et al. (2005) 


IC 3568 


0.54 


1.49E-16 


<0.0050 


<0.0039 


Liu et al. (2004) 


IC418 


1.43 


2.38E-15 


0.0046 


0.0005 


Henry, Kwitter & Bates (2000) 


IC 4406 


0.62 


2.37E-16 


<0.0022 


<0.0042 


Tsamis et al. (2003) 


IC4997 


0.34 


2.31E-16 


<0.002I 


<0.0046 


Hyung et al. (1994) 


IC5I17 


1.32 


8.29E-16 


0.0174 


0.0136 


Hyung et al. (2001) 


Ml-42 


0.15 


1.29E-16 


<0.0047 


<0.0018 


Liu et al .(2001b) 


M2-36 


0.76 


9.15E-17 


<0.0012 


<0.0084 


Tsamis et al. (2003) 


M4-18* 


2.9 


1.80E-16 


0.0491 


0.0111 


De Marco & Barlow (2001) 


Mz3 


0.83 


7.42E-15 


<0.0017 


<0.0003 


Zhang & Liu (2002) 


NGC 2346 


0.35 


1.15E-16 


<0.0087 


<0.0015 


Rola & Stasinska (1994); Peimbert & Torres-Peimbert (1987) 


NGC 2440* 


0.91 


4.49E-16 


<0.0043 


<0.0044 


Hyung & Aller (1998) 


NGC 39 18 


1.60 


9.59E-16 


<0.0033 


<0.0019 


Ercolano et al. (2003) 


NGC 40 


1.41 


1.38E-15 


0.0095 


0.0057 


Liu et al. (2004) 


NGC 5 189 


0.34 


1.71E-16 


<0.0018 


<0.0024 


Liu et al. (2001a) 


NGC 5315 


0.85 


1.12E-15 


0.0067 


0.0039 


Pottasch et al. (2002a) 


NGC 6153 


0.55 


1.20E-15 


<0.0022 


<0.0025 


Liu et al. (2000) 


NGC 6210 


0.26 


4.21E-16 


<0.0035 


<0.0059 


Liu et al. (2004) 


NGC 6302 


0.31 


8.08E-15 


0.0029 


0.0011 


Tsamis et al. (2003) 


NGC 6369 


2.00 


1.54E-15 


<0.0104 


<0.0055 


Zuckerman & Aller (1986);Aller & Keyes (1987) 


NGC 6445 


0.81 


3.38E-16 


<0.0026 


<0.0034 


Pottasch et al. (2002b) 


NGC 6537 


0.95 


1.50E-15 


0.0038 


0.0034 


Liu et al. (2001) 


NGC 6543 


0.44 


2.19E-15 


<0.0023 


<0.0003 


Wesson & Liu (2004) 


NGC 6572 


1.55 


1.55E-15 


<0.0076 


<0.0033 


Liu et al. (2004) 


NGC 6720 


0.62 


2.64E-16 


<0.0056 


<0.0008 


Liu et al. (2004) 


NGC 6741* 


1.35 


2.78E-16 


0.0070 


<0.0041 


Liu et al. (2004) 


NGC 6790 


0.81 


6.66E-16 


0.0129 


0.0015 


Liu et al. (2004) 


NGC 6826 


0.87 


1.05E-15 


<0.0077 


<0.0057 


Liu et al. (2004) 


NGC 6884* 


0.93 


3.03E-16 


0.0046 


0.0017 


Liu et al. (2004) 


NGC 6886* 


1.3 


2.50E-16 


<0.0048 


<0.0109 


Hyung etal. (1995) 


NGC 6891 


0.51 


1.62E-16 


<0.0158 


<0.0044 


Wesson et al. (2005) 


NGC 7009 


0.32 


1.13E-15 


<0.0011 


<0.0005 


Liu et al. (2004); Liu et al. (1995) 


NGC 7027 


2.29 


2.05E-15 


0.0171 


0.0106 


Middlemass (1990) 


NGC 7662 


1.23 


5.60E-16 


<0.0025 


<0.0016 


Liu et al. (2004) 


SwStl* 


0.93 


1.47E-15 


0.0083 


0.0018 


De Marco et al. (2001) 


Vy 2-2* 


0.26 


8.46E-16 


<0.0076 


<0.0062 


Wesson et al. (2004) 



error of the mean of 0.62±0.13 (13 PNe) and 0.40±0.06 (14 PNe), 
respectively. 

These empirical correlations between PAH band strength and 
nebular C/O abundance ratio could yield insights into the carbon 
dust condensation process because PAHs have been suggested to 
be the likely molecular precursors of C grains (Crawford, Tie- 
lens & AUamandola 1985; AUamandola, Tielens & Barker 1989). 
I(7.7)/I(IR) essentially measures the far-ultraviolet absorption cross 
sections of PAHs with respect to that of the larger C-grains, be- 
lieved to dominate the continuum emission in planetaries. One ob- 
vious interpretation of the observed correlations in Figs.|3|and|4|is 
that n(PAH)/n(carbon grains) increases with increasing C/O ratio. 
If, on the other hand, n(PAH)/;2(carbon grains) is largely invariant, 
one could interpret the correlations of PAH band strengths with 
C/O as due to growth of the grain component (higher C/O leading 



to bigger particles with smaller UV cross sections per unit volume, 
so that PAHs absorb more UV energy relative to the grains). 



7 ACKNOWLEDGEMENTS 

We thank the anonymous referee for comments that have helped 
improve the paper. MC thanks NASA for supporting the early por- 
tion of this work through grant NAG5-4884, and later through con- 
tract JPL96 1501, with Berkeley. It is a pleasure to thank the Physics 
& Astronomy Dept. at UCL for hosting a number of short-term vis- 
its that have enabled the completion of this effort. 



PAHs and C/0 ratios in planetaries 9 



Table 3. Detected PAH 6.2-//m ratios, I(6.2)/I(IR), for PNe 



PN 


I(6.2)/I(IR) 




0.0073 


Cnl-5 


0.0041 


CPD-56°8032 


0.0105 


He 2-1 13 


0.0116 


IC418 


0.0021 


IC5117 


0.0140 


M4-18 


0.0051 


NGC5315 


0.0037 


NGC6302 


0.0008 


NGC6790 


0.0089 


NGC6826 


0.0024 


NGC 6886 


0.0079 


NGC 7027 


0.0039 


SwStl 


0.0050 
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